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Abstract 

Supersymmetric models which fulfill the conditions of electroweak naturalness generally 
contain light higgsinos with mass not too far from mh — 125 GeV, while other sparticles 
can be much heavier. In i?-parity conserving models, the lightest neutralino is then 
a higgsino-like WIMP (albeit with non-negligible gaugino components), with thermal 
relic density well below measured values. This leaves room for axions to function as 
co-dark matter particles. The local WIMP abundance is then expected to be below 
standard estimates, and direct and indirect detection rates must be accordingly rescaled. 
We calculate rescaled direct and indirect higgsino-like WIMP detection rates in SUSY 
models that fulfill the electroweak naturalness condition. In spite of the rescaling, we find 
that ton-scale noble liquid detectors can probe the entire higgsino-like WIMP parameter 
space, so that these experiments should either discover WIMPs or exclude the concept 
of electroweak naturalness in i?-parity conserving natural SUSY models. Prospects for 
spin-dependent or indirect detection are more limited due in part to the rescaling effect. 
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1 Introduction 



The recent discovery of a Higgs-like resonance by Atlas and CMS collaborations [TJ [2] at ~ 
125 GeV seemingly adds another supporting pillar to the theory of weak scale supersymmetry, 
since in the Minimal Supersymmetric Standard Model (MSSM) [3] we expect rrih — 115 — 135 
GeVjl]. Aside from solving the gauge hierarchy problem, previous supporting pillars arising 
from data include 1. the measured values of gauge couplings allow for gauge coupling unification 
at a scale Q ~ 2 x 10 16 GeV within the MSSM, 2. the large value of the top quark mass is 
precisely what is needed to drive electroweak symmetry breaking via radiative corrections and 3. 
SUSY is replete with several possible cold dark matter (CDM) candidates (neutralino/WIMP, 
gravitino, axino) in the form of the lightest SUSY particle (LSP). 

Since a variety of SUSY model parameters enter into the scalar potential of the theory, and 
thus contribute to the W, Z and Higgs masses, it is widely expected that superpartners should 
exist at or around the weak scale. This mantra has been repeated in numerous talks and papers 
over the past decades, so we refer to it here as the story of SUSY electroweak naturalness^. 
Indeed, the concept of naturalness may dictate to some degree when it is time to give up on 
weak scale supersymmetry should no signal be ultimately found[5]. 

While the existence of a Higgs-like scalar at ~ 125 GeV is a boon for SUSY models, on the 
contrary, no signal for superpartners has yet emerged at LHC. In models such as the popular 
mSUGRA/CMSSM[6], the Atlas and CMS collaborations [8] now require m g ~ 1.4 TeV for 
trig ~ rrig and m g ~ 1 TeV for rrig 3> m g . Already this fact has led some astute physicists to 
give up on weak scale SUSYj9], or to at least concede that weak scale SUSY is finetuned. Thus, 
the recent LHC limits on sparticle masses seemingly exacerbate what is known as the Little 
Hierarchy Problem (LHP): why is there such a disparity between the sparticle mass scale and 
the electroweak scale? 

Before jumping to conclusions, it pays to scrutinize electroweak naturalness more closely. 
We can be more precise if we re-phrase the LHP in the following terms: how is it that the 
Z-boson mass can exist at just 91.2 GeV while gluino and squark masses are at, or even well 
beyond, the TeV scale? The answer proposed in Ref 's JTUJ [TlJ [12] is that all the individual 
weak scale contributions feeding mass into rriz should be not too far from rriz- The value of 
rriz in the MSSM is given by 



m% _ (m^ + S^)-(m^ + S^)tan 2 /3 2 2 
2 " (tan 2 /3-l) " ^ 



2 2 2 I -i \ 



where the latter approximate equality obtains for ratio-of-Higgs vevs tan/3 = v u /vd ~ 3 and 
where the and terms represent the sum of various radiative correct ions [12]. To be 
quantitative, a fmetuning measure 



±EW = rnaxi 



Ci/(m 2 z /2) (2) 



may be defined, where Ci represents any of the terms on the right-hand-side of Eq. [T] {e.g. 
Ch u = — m H u tan 2 /3/(tan 2 /3 — 1) and = — /i 2 ). The fmetuning measure A EW enjoys several 
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advantages as discussed in Ref. [12]: it is 1. model-independent (within the MSSM) in that 
any model giving rise to the same weak scale SUSY spectrum will have the same measure of 
finetuning, 2. conservative, 3. unambiguous, 4. predictive, 5. falsifiable and 6. simple to 
encode. 

In the 2-parameter non-universal Higgs model (NUHM2)[13], with parameter space given 

by 

m , mi/2, Aq, tan/3, fj,, m A , (3) 

scans have found that values of A EW as low as ~ 10 (corresponding to A EW ~ 10% electroweak 
finetuning) could be found[10l [12]. To achieve low A E w, one needs a) \fj,\ ~ mz ~ 100 — 200 
GeV, b) mn u ~ (1 — 2)mo so that m% is driven radiatively (via the large top-quark Yukawa 
coupling) to small but negative values at the weak scale and c) moderate values of m 1( / 2 ~ 
0.3 — 1 TeV. Scalar masses may exist in the range 1-10 TeV (even higher if one allows for split 
generations) although large trilinear soft breaking terms A ~ 1.5m are required. The large 
A-terms yield large mixing in the top-squark sector which simultaneously softens the radiative 
corrections S^(ti) and £^(£2), while lifting the value of rrih into the 125 GeV range [TO] . While 
the NUHM2 model admits values of A EW even below 10, the mSUGRA/CMSSM model is 
considerably more tuned, where a minimum A E w ~ 100 has been found, although typically 
A E w f° r mSUGRA is more like 10 3 — 10 4 . Models with low A E w ~ 30 have been dubbed 
radiative natural SUSY (RNS) since the low electroweak finetuning is radiatively driven. 
The sparticle spectra found for RNS models with A E w ~ 30 is characterized by: 

• light higgsino-like and Z 12 with masses ~ 100 — 300 GeV, 

• gluinos with mass m g ~ 1 — 4 TeV, 

• top squarks with mi ~ 1 — 2 TeV and m^ 2 ~ 2 — 5 TeV, 

• first / second generation squarks and sleptons with masses m~ ~ t ~ 1 — 8 TeV. The m~ ~ t range 
can be pushed up to 20-30 TeV if non-universality of generations with mo(l, 2) > m (3) 
is allowed. 

The RNS model with the above spectra yields branching fractions BF{b — > 57) and BF(B S — > 
/i + /i~) in accord with measured values, unlike many models with lighter top squarks [T4]. 

As far as testability goes, RNS models should yield observable signals from gluino pair 
production at LHC14 with ~ 300 fb _1 for m- g ~ 1.6 TeV [15], The LHC14 reach for SS diboson 
production from pp — > W^Z^ — > W ± W ± + $t gives a reach in terms of m g of m g ~ 1.8 TeV 
for the same integrated luminosity [16]. Since m g can range up to ~ 4 TeV for the RNS models, 
then LHC14 will not be able to definitively discover /exclude RNS. 

Alternatively, the hallmark distinction of RNS models is the presence of four light higgsinos 
with masses ~ 100 — 300 GeV. A linear e + e~ collider operating at y/s ~ 2\fi\ can definitively 
discover /exclude RNS models, e.g. ^J~s ~ 0.5 (1) TeV would access all models with A E w ~ 25 
(70). Since the timescale for operation of a TeV-scale ILC is of order 10-20 years, here we 
address instead the possibility of much earlier discovery of the higgsino-like WIMPs expected 
from RNS models. 
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2 Direct and indirect higgsino detection 



In this paper, we generate RNS models from the Isasugra spectrum generator [T7] using the 
same NUHM2 parameter space scan as utilized in Ref. [12]. We will require that the model 
generated from each parameter set obeys the LEP2 bound that > 103.5 GeV, and will focus 
on models with Aew < 50 (100) corresponding to better than 2% (1%) electroweak fmetuning. 
We calculate the "standard thermal neutralino abundance" Q~ d h 2 using the IsaReDfl8l relic 

n 

density subroutine^ We will accept only models with Q~ d h 2 < 0.12 for reasons to be made 
clear shortly. The relic abundance from RNS models is shown in Fig. [TJ The red crosses have 
Aew < 50 whilst blue dots have Aew < 100. Green points we will find later are already 
excluded by direct /indirect WIMP search limits. From the figure, we see a high density band 
extending from Q^h 2 ~ 0.004 for m 5 ~ 100 GeV to Vt^h 2 ~ 0.02 for m= ~ 300 GeV, i.e. 
there is typically a standard underabundance of higgsino dark matter compared to measurement 
from WMAP9|19j by a factor ranging from 3-25. There is some spread in these values above 
and below the main band from cases where /i is quite large and mi/ 2 is small so that one 
has instead a mixed higgsino-bino LSP state. The bulk of points above the band are already 
excluded as we shall see. Thus, the mainly higgsinolike neutralino by itself does not make a 
good CDM candidate. Additional new physics is needed to match the measured dark matter 
density. 

One compelling way to address the dark matter deficiency is by invoking the Peccei-Quinn- 
Weinberg-Wilczek solution to the strong CP problemJ^O] via introduction of a PQ symmetry 
and concommitant axions. Since we are working in SUSY models, the axion will be accompanied 
by i?-parity-even spin-0 saxions s and .R-parity-odd spin-| axinos a[21]. In gravity-mediation 
models (as assumed here), the saxion and axino are expected to have masses m s ~ m s ~ 
m 3/2 ~ 5 — 20 TeV[22j. the same mass scale as matter scalars m in the theory. In this case, 
the dark matter would consist of both axions and higgsinos acting as co-dark matter particles. 

The relic abundance of mixed axion-neutralino CDM has been addressed in Ref 's [231 EH 
[25| [26] . In [26] , it was found that SUSY models with a standard overabundance of dark matter 
are still excluded in the PQMSSM by a combination of dark matter overabundance constraints, 
BBN constraints and dark radiation constraints. However, SUSY models with a standard 
underabundance of neutralinos are still allowed over large ranges of PQMSSM parameters. For 
models with a standard underabundance, then thermal production and decay of axinos in the 
early universe augments the neutralino abundance, sometimes by too much, other times not 
enough: the former case would be excluded whilst in the latter case, the remaining abundance 
is made up of relic axions produced through the usual vacuum misalignment mechanism. In 
addition, for large values of PQ breaking scale f a ~ 10 13 — 10 16 GeV, then saxions can be 
produced at large rates via coherent oscillations. The saxions can augment the neutralino 
abundance by decaying to SUSY particles (e.g. s — > gg), or they can dilute all relics if they 
decay after freezeout into SM particles. Saxions may also decay into axion pairs s — > aa leading 
to production of dark radiation [27]. From the PQMSSM scans in Ref. [26], it is found that 
in the case where s — > aa branching fraction is at all substantial, then entropy dilution is 

2 The standard thermal abundance refers to the neutralino density derived from assuming neutralinos in 
thermal equilibrium at high temperature followed by freeze-out due to the Universe's expansion. 
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■ a h =0.11 




Figure 1: Plot of standard thermal neutralino abundance Q~ d h 2 versus m(higgsino) from a 
scan over NUHM2 parameter space with Aew < 50 (red crosses) and Aew < 100 (blue dots). 
Green points are excluded by current direct /indirect WIMP search experiments. We also show 
the central value of Q C DMh 2 from WMAP9. 
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always accompanied by a violation of either BBN or dark radiation constraints (parametrized 
by AN e fj ~ 1.6 where AN e ff is the effective number of additional neutrinos beyond the SM 
value). Thus, the scans over PQMSSM parameter space in Ref. [25] find that the standard 
underabundance can be augmented by any factor leading to Q~ d h 2 < Q^Jt 2 < 0.11, but not 
diminished without violating BBN or DR constraints. Models with too much CDM production, 
or models which violate BBN or dark radiation constraints would be excluded. The upshot is 
that in RNS models, for any particular parameter set, we expect the relic higgsino abundance to 
lie somewhere between the standard value Q~ d h 2 (which would correspond to axion domination) 
up to the measured value 0.11, in which case CDM would be higgsino-dominated. The question 
then arises: what are the prospects for direct /indirect detection of relic higgsinos in WIMP 
detection experiments? 

In Fig. [21 we show the spin-independent higgsino-proton scattering rate in pb as calculated 
using IsaReS[2S]. The result is rescaled by a factor £ = Q~ d h 2 /0.11 to account for the fact that 

the local relic abundance might be far less than the usually assumed value pi oca i — 0.3 GeV/cm 3 , 
as suggested long ago by Bottino et a/.[2H] (the remainder would be composed of axions). The 
higgsino-like WIMP in our case scatters from quarks and gluons mainly via h exchange. The 
Z\ — Zi — h coupling involves a product of both higgsino and gaugino components. In the 
case of RNS models, the Z\ is mainly higgsino-like, but since m\n is bounded from above by 
naturalness, the Z\ contains enough gaugino component that the coupling is never small: in 
the notation of Ref. [3] 

C 3 -X^Z^h (4) 

where 

Vh 1 / (1) ■ (1) \ ( (1) / ,r\ 

X n = -- [yi 'sina-Dj' cos a ) [gv$ - gv\'j , (5) 

and where and v^p are the higgsino components and and vjp are the gaugino compo- 
nents of the lightest neutralino, a is the Higgs mixing angle and g and g' are SU(2)l and U(l)y 
gauge couplings. Thus, for SUSY models with low Aew ~ 50 — 100, the SI direct detection 
cross section is also bounded from below, even including the rescaling factor £. 

From the Figure, we see that the current reach from 225 live-days of Xe-100 running|30j 
already bites into a significant spread of parameter points. The excluded points are colored 
green. The projected reach of the LUX 300 kg detector[31] is also shown by the black-dashed 
contour, which should explore roughly half the allowed RNS points. The reach of SuperCDMS 
150 kg detcctor[32j is shown as the purple-dashed contour. The projected reach of Xe- 1-ton, 
a ton scale liquid Xenon detector, is also shown[33j. Our main result is this: the projected 
Xe-l-ton detector, or other comparable noble liquid detectors, can make a complete exploration 
of the RNS parameter space. Since deployment of the Xe-l-ton detector is imminent, it seems 
direct WIMP search experiments may either verify or exclude RNS models in the near future, 
thus bringing the story of electroweak naturalness to a conclusion! 

While the above result is indeed compelling, it is not a theorem, and is subject to several 
reasonable assumptions. These include the assumption of i?-parity conservation with a higgsino- 
like co-DM particle, along with the assumption of non-negligible s — > aa decay rate. If the 
s -y od decay rate is somehow forbidden or highly suppressed, then s — >• gg (gluons) may be 
dominant, in which case substantial entropy dilution can still occur [24"l [25] . 
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m(higgsino) (GeV) 



Figure 2: Plot of rescaled higgsino-like WIMP spin- independent direct detection rate £a SI (Zip) 
versus m(higgsino) from a scan over NUHM2 parameter space with Aew < 50 (red crosses) 
and Aew < 100 (blue dots). Green points are excluded by current direct /indirect WIMP 
search experiments. We also show the current reach from Xe-100 experiment, and projected 
reaches of LUX, SuperCDMS 150 kg and Xe-1 ton. 
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Figure 3: Plot of rescaled spin-dependent higgsino- like WIMP detection rate ^a SD (Z 1 p) versus 
m(higgsino) from a scan over NUHM2 parameter space with Aew < 50 (red crosses) and 
A-ew < 100 (blue dots). Green points are excluded by current direct /indirect WIMP search 
experiments. We also show current reach from the COUPP detector. 

An alternative method for rescuing theories with an underabundance of WIMPs is to hy- 
pothesize the existence of some scalar field (such as a modulus field from string theory) which 
can be produced via coherent oscillations, and which can engage in late decays injecting ei- 
ther more WIMPs (thus increasing the WIMP abundance) or entropy (thus decreasing the 
abundance [31]). The resulting abundance just depends on two parameters: the scalar field 
decay temperature and the branching fraction into SUSY particles [35], [361 [33, [38]. In our case, 
with a higgsino-like WIMP underabundance, the decaying scalar field would increase the hig- 
gsino abundance to its measured value, and the rescaling factor £ = 1. For this possibility, the 
results of Ref. [12] would apply. 

In Fig. (3[ we show the rescaled spin- dependent higgsino-proton scattering cross section 
£<t (Zip). Here we show recent limits from the COUPP [39] detector. Current limits are still 
about an order of magnitude away from reaching the predicted rates from RNS models. 

To compare against the current reach of IceCubejlO], we show in Fig. S] the value of 
a SD (Zip), with no rescaling factor. Here, the IceCube rates should not be rescaled since the 
IceCube detection depends on whether the Sun has equilibrated its core abundance between 
capture rate and annihilation rate[JT]. Typically for the Sun, equilibration is reached for almost 
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Figure 4: Plot of (non-rescaled) spin-dependent higgsino- like WIMP detection rate a SD {Z 1 p) 
versus m(higgsino) from a scan over NUHM2 parameter space with Aew < 50 (red crosses) 
and A EW < 100 (blue dots). Green points are excluded by current direct /indirect WIMP 
search experiments. We also show current reach from IceCube. 

all of SUSY parameter space [32] ■ The IceCube limits have entered the RNS parameter space 
and excluded the largest values of a SD (Zip). 

In Fig. we show the rescaled thermally-averaged neutralino annihilation cross section 
times relative velocity in the limit as v — > 0: ^ 2 (crv ) |„_>o- This quantity enters into the rate 
expected from WIMP halo annihilations into 7, e + , p or d. The rescaling appears as £ 2 since 
limits depend on the square of the local WIMP abundance [33]. Anomalies in the positron and 
7 spectra have been reported, although the former may be attributed to pulsars[44], while the 
latter 130 GeV gamma line may be instrumental. Soon to be released results from AMS-02 
should clarify the situation. On the plot, we show the limit derived from the Fermi LAT gamma 
ray observatory[45j for WIMP annihilations into WW. These limits have not yet reached the 
RNS parameter space due in part to the squared rescaling factor. 
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Figure 5: Plot of rescaled ^ 2 (cv)\ v ^q versus m(higgsino) from a scan over NUHM2 parameter 
space with A EW < 50 (red crosses) and A EW < 100 (blue dots). Green points are excluded 
by current direct /indirect WIMP search experiments. We also show current reach from Fermi 
LAT, Ref. H5J. 
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3 Conclusions: 



In conclusion, we have found that SUSY models can elude the Little Hierarchy Problem in 
the guise of radiative natural SUSY models which feature a mainly higgsino-like neutralino 
that may act as a co-dark-matter particle along with the axion. While LHC can explore a 
portion of RNS parameter space, an ILC can probe it entirely, although such a machine may 
be well over a decade in the future. However, current WIMP direct detection experiments 
are biting into the meat of RNS parameter space, even if we take into account rescaling of 
the local abundance due to the fact that higgsinos may make up only a portion of the dark 
matter. LUX and ultimately SuperCDMS will probe further. The soon-to-be-deployed Xe- 
1-ton noble liquid detector should be able to completely explore the entire RNS parameter 
space (as shown in Fig. [5]), thus either discovering a higgsino-like WIMP or rejecting the story 
of SUSY electroweak naturalness. Complementary signals from spin-dependent and indirect 
WIMP detection channels are less likely since these are usually suppressed by the reduced local 
WIMP density which is expected from models of mixed dark matter. 

Acknowledgments 

We thank K. J. Bae, Y. Gershtein, P. Huang, A. Lessa, A. Mustafayev and X. Tata for dis- 
cussions. This work was supported in part by the US Department of Energy, Office of High 
Energy Physics. 

References 

[1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716 (2012) 1. 

[2] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716 (2012) 30. 

[3] H. Baer and X. Tata, Weak Scale Supersymmetry: From Superfields to Scattering Events, (Cam- 
bridge University Press, 2006). 

[4] M. S. Carena and H. E. Haber, Prog. Part. Nucl. Phys.50 (2003) 63 [hep-ph/0208209]. 

[5] For a recent review, see e.g. J. L. Feng, arXiv: 1302.6587 [hep-ph]. 

[6] A. Chamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 49 (1982) 970; R. Barbieri, S. Ferrara 
and C. Savoy, Phys. Lett. B 119 (1982) 343; N. Ohta, Prog. Theor. Phys. 70 (1983) 542; L. Hall, 
J. Lykken and S. Weinberg, Phys. Rev. D 27 (1983) 2359; for a recent review, see R. Arnowitt, 
A. H. Chamseddine and P. Nath, Int. J. Mod. Phys. A 27 (2012) 1230028. 

[7] G. Aad et al. (ATLAS collaboration), Phys. Lett. B 710 (2012) 67. 

[8] S. Chatrchyan et al. (CMS collaboration), Phys. Rev. Lett. 107 (2011) 221804. 

[9] See e.g. M. Shifman, Mod. Phys. Lett. A 27 (2012) 1230043. 

[10] H. Baer, V. Barger, P. Huang, A. Mustafayev and X. Tata, Phys. Rev. Lett. 109 (2012) 161802. 



10 



[11] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev and X. Tata. larXiv: 1210.30191 

[12] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev and X. Tata, arXiv:1212.2655l [hep- 
ph]. 

[13] J. Ellis, K. Olive and Y. Santoso, Phys. Lett. B 539 (2002) 107; J. Ellis, T. Falk, K. Olive and 
Y. Santoso, Nucl. Phys. B 652 (2003) 259; H. Baer, A. Mustafayev, S. Profumo, A. Belyaev and 
X. Tata, J. High Energy Phys. 0507 (2005) 065. 

[14] H. Baer, V. Barger, P. Huang and X. Tata, J. High Energy Phys. 1205 (2012) 109. 

[15] H. Baer, V. Barger, A. Lessa and X. Tata, Phys. Rev. D 86 (2012) 117701. 

[16] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev, W. Sreethawong and X. Tata, 
larXiv:1302.58TBl [hep-ph] . 

[17] H. Baer, C. H. Chen, R. Munroe, F. Paige and X. Tata, Phys. Rev. D 51 (1995) 1046; H. Baer, 
J. Ferrandis, S. Kraml and W. Porod, Phys. Rev. D 73 (2006) 015010. 

[18] IsaReD, see H. Baer, C. Balazs and A. Belyaev, J. High Energy Phys. 0203 (2002) 042. 

[19] G. Hinshaw, D. Larson, E. Komatsu, D. N. Spergel, C. L. Bennett, J. Dunkley, M. R. Nolta and 
M. Halpern et al. larXiv:1212.5226l . 

[20] R. Peccei and H. Quinn, Phys. Rev. Lett. 38 (1977) 1440 and Phys. Rev. D 16 (1977) 1791; S. 
Weinberg, Phys. Rev. Lett. 40 (1978) 223; F. Wilczek, Phys. Rev. Lett. 40 (1978) 279. 

[21] H. P. Nilles and S. Raby, Nucl. Phys. B 198 (1982) 102; J. E. Kim, Phys. Lett. B 136 (1984) 
378; J. E. Kim and H. P. Nilles, Phys. Lett. B 138 (1984) 150; for a review, see e.g. F. D. Steffen, 
Eur. Phys. J. C 59 (2009) 557. 

[22] P. Moxhay and K. Yamamoto, Phys. Lett. B 151 (1985) 363; E. Chun and A. Lukas, Phys. Lett. 
B 357 (1995) 43; J. E. Kim, M.-S. Seo, Nucl.Phys.B864 (2012) 296 ; C. Cheung, G. Elor and 
L. J. Hall, Phys. Rev. D 85 (2012) 015008. 

[23] K-Y. Choi, J. E. Kim, H. M. Lee and O. Seto, Phys. Rev. D 77 (2008) 123501. 

[24] H. Baer, A. Lessa, S. Rajagopalan and W. Sreethawong, JCAP1106 (2011) 031. 

[25] H. Baer, A. Lessa and W. Sreethawong, JCAP1201 (2012) 036. 

[26] K. J. Bae, H. Baer and A. Lessa. llarXiv:1301.7428l [hep-ph], 

[27] See e.g. V. Barger, J. P. Kneller, H. -S. Lee, D. Marfatia and G. Steigman, Phys. Lett. B 
566 (2003) 8; K. Ichikawa, M. Kawasaki, K. Nakayama, M. Senami and F. Takahashi, JCAP 
0705 (2007) 008; J. Hasenkamp, Phys. Lett. B 707 (2012) 121; J. Hasenkamp and J. Kersten, 
larXiv:1212.4l"60 ; D. Hooper, F. Queiroz and N. Gnedin, Phys. Rev. D 85 (2012) 063513; E. Di 
Valentino, S. Galli, M. Lattanzi, A. Melchiorri, P. Natoli, L. Pagano and N. Said, arXiv: 1301. 7343 
[astro-ph.CO]; M. Cicoli, J. P. Conlon and F. Quevedo, larXiv: 1208.35621 [hep-ph]; P. Graf and 
F. D. Steffen. larXiv:1208.295ll and larXiv:1302.2143l [hep-ph]. 

[28] H. Baer, C. Balazs, A. Belyaev and J. O'Farrill, JCAP 0309 (2003) 007. 



11 



[29 
[30 
[31 

[32 
[33 
[34 
[35 

[36 
[37 
[38 

[39 
[40 
[41 
[42 
[43 
[44 



A. Bottino, F. Donato, N. Fornengo and S. Scopel, Phys. Rev. D 63 (2001) 125003. 

E. Aprile et al. [XENON100 Collaboration], Phys. Rev. Lett. 109 (2012) 181301. 

D. S. Akerib et al. [LUX Collaboration], Nuclear Inst. Methods in Physics Research A704 (2013) 
111-126 [arXiv:1211.3788l [Dhvsics.ins-det]]. 



P. L. Brink et al. [CDMS-II Collaboration], eConf C 041213 (2004) 2529 |astro-ph/0503583| . 
E. Aprile [XENON1T Collaboration]. llarXiv: 1206.62881 
T. Moroi and L. Randall, Nucl. Phys. B 570 (2000) 455. 

G. Gelmini and P. Gondolo, Phys. Rev. D 74 (2006) 023510; G. Gelmini, P. Gondolo, A. Sol- 
datenko and C. Yaguna, Phys. Rev. D 74 (2006) 083514. 

M. Endo and F. Takahashi, Phys. Rev. D 74 (2006) 063502. 

B. Acharya, G. Kane, S. Watson and P. Kumar, Phys. Rev. D 80 (2009) 083529. 

R. Allahverdi, B. Dutta and K. Sinha, Phys. Rev. D 86 (2012) 095016; R. Allahverdi, B. Dutta 
and K. Sinha, larXiv: 1212.69481 [hep-ph]. 

E. Behnke et al. [COUPP Collaboration], Phys. Rev. D 86 (2012) 052001. 

R. Abbasi et al. (IceCube collaboration), Phys. Rev. D 85 (2012) 042002. 

G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267 (1996) 195. 

V. Niro, A. Bottino, N. Fornengo and S. Scopel, Phys. Rev. D 80 (2009) 095019. 

A. Bottino, F. Donato, N. Fornengo and P. Salati, Phys. Rev. D 72 (2005) 083518. 

V. Barger, Y. Gao, W. Y. Keung, D. Marfatia and G. Shaughnessy, Phys. Lett. B 678 (2009) 
283 ; S. Profumo, Central Eur. J. Phys. 10 (2011) 1. 



[45] M. Ackermann et al. (Fermi Collaboration), Phys. Rev. Lett. 107 (2011) 241302; A. Geringer- 
Sameth and S. M. Koushiappas, Phys. Rev. Lett. 107 (2011) 241303. 



12 



